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We investigate how well future large-scale radio surveys could measure different shapes of primordial non- 
Gaussianity; in particular we focus on angle-dependent non-Gaussianity arising from primordial anisotropic 
sources, whose bispectrum has an angle dependence between the three wavevectors that is characterized by 
Legendre polynomials Vl and expansion coefficients cl . We provide forecasts for measurements of galaxy 
power spectrum, finding that Large-Scale Structure (LSS) data could allow measurements of primordial non- 
Gaussianity competitive or improving upon current constraints set by CMB experiments, for all the shapes 
considered. We argue that the best constraints will come from the possibility to assign redshift information 
to radio galaxy surveys, and investigate a few possible scenarios for the EMU and SKA surveys. A realistic 
(futuristic) modeling could provide constraints of /nl ~ 1(0.5) for the local shape, /nl of O(10)(O(l)) for 
the orthogonal, equilateral and folded shapes, and cl= i ~ 80(2), cl= 2 ~ 400(10) for angle-dependent non- 
Gaussianity. The more futuristic forecasts show the potential of LSS analyses to considerably improve current 
constraints on non-Gaussianity, and so on models of the primordial Universe. Finally, we find the minimum 
requirements that would be needed to reach ct(cl= 1 ) = 10, which can be considered as a typical (lower) value 
predicted by some (inflationary) models. 


I. INTRODUCTION 


Investigating models describing the primordial universe is a priority for current and planned cosmological experiments; recent 
results from measurements of the CMB by e.g. the Planck satellite OJ0 and of the Large-Scale Structure (LSS) of the Universe, 
by e.g. the SDSS EM) already set some constraints on such models, but more precision is needed in order to further discriminate 
between them. The window for investigating the inflationary period is given by the statistics of primordial curvature fluctuations, 
and in particular the (deviations from) gaussianity of the probability distribution function of primordial cosmological perturba¬ 
tions. Imprints of primordial non-Gaussianity in the CMB anisotropies and in the LSS are usually described in terms of shapes 
(e.g. EHHl)- Recently, along with “standard” shapes like the local, equilateral, orthogonal and folded templates, shapes with 
specific angle dependence have also been investigated (e.g., EMI). 

The angle dependence of primordial non-Gaussianity can be generated by the existence of non-Gaussian anisotropic sources. 
We may parameterize the angular dependence in the bispectrum of curvature perturbations using Legendre polynomials 
Pl(x) Go): 


B^(k 1 ,k 2 ,k 3 ) 


cl 


Vl (ki • k 2 )^c( fc i)- P c( fc 2 ) + 2 perm , 


( 1 ) 


where Co is related to the local-type non-Gaussianity parameter as Co = (6/5)/^ • Non-Gaussian anisotropic sources actually 
create nonzero cl for L > 1. For example, in a model where the inflaton field c\> is coupled to a U(l) gauge field via a 
f{(j))F 2 interaction (where F^ v = d^A u — d u A^ is the vector field strength), non-vanishing C 2 arise as well as Co, whose 
magnitudes are related to each other as C 2 = co/2 [18] 2H]. 1 The case of C 2 Co can be realized in the so-called solid 
inflation models |[T9l[2Tl[24ll25l . which are based on a specific internal symmetry obeyed by the inflaton fields. In the presence 
of large-scale non-helical and helical magnetic fields in the radiation-dominated era, not only Co and C 2 but also c\ can be 
created dunum Recently a model with a /(</>) (F 2 + FF) coupling has been proposed as the first clear example of an 
inflationary model where c\ is generated, realizing cq : c\ : C 2 = 2 : — 3 : 1 l26l . 


1 Notice in the f(4>)F 2 where the vector field has a non-vanishing vev a statistical anisotropic bispectrum is actually generated and, after an angle-average, it 
takes the form |TJ. For studies of bispectra that break statistical isotropy see I22II23I . 
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Therefore, c 1 and C 2 can become key observables to probe these models. So far, observational constraints have been investi¬ 
gated via the CMB bispectrum analysis. In an all-sky ideal measurement of the temperature anisotropies up to £ = 2000, where 
instrumental noise is completely negligible, expected lcr error bars on c\ and C 2 reach 61 and 13, respectively (20) . Realistic 
constraints have been obtained from the Planck temperature data, giving c\ — 118 =b 103 and C 2 = — 5 ± 26 (68% CL). 2 

As well as in the CMB I4l[5lfl0l [271 . signatures of primordial non-Gaussianity can be also measured in the LSS of the Universe 
via the scale dependence of the halo bias (see e.g. l28ti30l ). In this paper, we examine how powerful can LSS surveys be in 
measuring the angle-dependent non-Gaussianity parameters c\ and C 2 which could reveal about specific details of the physics of 
inflation, such as the presence of vector fields. 3 

Galaxy clustering can and has been used to set constraints on e.g. dark energy parameters 1(331 . models of gravity (34), neutrino 
mass |35|[36), dark matter models (371(38) . the growth of structures (39ll40l . and here we are in particular interested in their 
use to measure primordial non-Gaussianity. Non-Gaussianity parameters can be measured via the LSS by using the Integrated 
Sachs-Wolfe effect 14111421 and the galaxy power spectrum (see e.g. E3 EH EH0309; for some forecasts on measurements 
on non-Gaussianity with 21-cm surveys, see e.g. EoHU. For a recent review on testing inflation with galaxy surveys, see in 
particular E3I3I). 

Recent measurements from LSS include EH ED; with the current generation of experiments, these constraints are not 
competitive with CMB ones, but future galaxy surveys will provide extremely precise measurements of galaxy clustering over a 
wide area of the sky and a large range in redshifts (see [53- 59). For this reason, planned experiments such as SPHEREx [ 601 will 
be explicitly focused on investigating the primordial universe (see also Eim for a general overview of expected constraints). 

The paper is organized as follows. In Section [II] we study the effects of angle-dependent non-Gaussianity on the halo bias 
and compare them to effects due to other shapes. In Section [HI] we present the survey specifications and assumptions and in 
Section [IV] we describe the observable used to compute our forecasts. In Section [V] we present the methodology adopted for 
obtaining our results and we present them, with a focus on future constraints on vector fields models in Section |VI| Finally, in 
Section IVTlI we draw our conclusions. 


II. HALO BIAS IN THE PRESENCE OF ANGLE-DEPENDENT NON-GAUSSIANITY 


It is well-known that non-vanishing non-Gaussianity can create nontrivial scale dependence in the bias parameter associating 
halos with matter fluctuations, in addition to a usual scale-invariant contribution in the Gaussian case bo- In the presence of 
nonzero primordial non-Gaussianity, one may write the total bias as: 


6ng0, k ) = frcO) + k ) • 


( 2 ) 


There are diverse approaches to estimating the scale-dependent part A b for given primordial non-Gaussianities (see e.g. [29 , 
[30ll4lll43ll44ll6ll ). In this work, we estimate A b by means of the iPT formalism (61), which is applicable to the analysis of 
almost any type of primordial non-Gaussianity. The corresponding formula reads: 


A b(z, k ) 


251 


A 2 (M)l(k,M)+A 1 (M) 


dl{k, M) 
d In (Jm 


( 3 ) 


where S c = 1.686 is the critical overdensity and gm denotes the density variance per the mass of halo, M 


^PmoR 3 ' 


= j ^ W\kR)P h {k ), 


(4) 


p m o being the mean matter density at the present epoch and Ph{k) the linear matter power spectrum. 4 As a window function, we 
adopt a top-hat form: W(x) = 3jl . The coefficients Ai(M) are determined by the mass function of halos. In the numerical 
works of this paper, we use fitting formulae for Ai(M) based on the MICE mass function (64) . whose explicit forms are derived 
in (61). 


2 These values correspond to the Planck 2015 limits 0: f^ 1 = —49 =b 43 and /= 0.5 d= 2.7, where the /^ L parameters in the Planck papers |[4|[5] 
are related to cl inEq. |lj as ci = — (12/5 )/j5(l 1 and C 2 = — (48/5) 2 . 

3 This paper focuses on the isotropic angle-averaged bispectrum effects on the scale-dependent bias due to primordial vector fields. See ED for the analysis 
including the statistically-anisotropic effects. For the analysis of the effects sourced from large-scale vector fields existing in the radiation dominated era, 
see f32l l. 

4 In this paper we only take into account the contributions from the bispectrum. See (22122) f° r discussions on the higher-order effects. 
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The function X is given by a convolution of the linear matter bispectrum B l (filtered with the window function): 

1 f k 7 

m M) = J j2^ w{k ' R)w ^ k ~ k '\ R ^ k ’ fc, > i k - k 'D • 


(5) 


In linear perturbation theory, matter fluctuations become a simple product of primordial curvature perturbation and the linear 
matter transfer function, i.e. S( k, z) = M^(k, z)(k, and hence their power spectrum and bispectrum at given z are written 
respectively as: 


P h (k) = M 2 c (k,z)P < (k), B h (k !,k 2 ,k 3 ) 


" 3 

II M^{k n , z) 

_n=1 


B c (k 1 ,k 2 ,k 3 ). 


In the standard ACDM cosmology: 


M ( (k,z) 


2 k 2 T(k) 


( 6 ) 


(7) 


a derivation can be found in Il32ll65ll66l . From these relations, one can deduce that the k dependence of A b varies according to 
the shape of the primordial bispectrum B £, via the variation of the function X. In the next section, by analyzing matter power 
spectra and the bias parameter, we will explore the detectability of angle-dependent non-Gaussianity, i.e., c\ and C 2 of Eq. ([!]), 
as well as the detectability of the usual four types of non-Gaussianity, i.e., /§l and /^. The data of A b used in such 

analysis are obtained through numerical computations of the iPT formalism as in Eq. Q. 

Our numerical results for A b(z, k) for each primordial non-Gaussianity model at z = 1 are shown in Figure [I] It is visually 
apparent there that the tilt of k varies based on the non-Gaussianity shape. Some simple scaling arguments can help in under¬ 
standing the main features of the results of Figure [T] It is a known fact that the convolution in the X function is determined, in the 
large-scale limit, by the behavior of B^ in the squeezed configurations, k k' ~ |k — k'|. In particular for the squeezed-limit 
of the the angle-dependent primordial bispectrum we find: 

B% =0 (k,k', |k-k'|) -► 2c 0 P c (k)P c (k') <xB l £ c (k,k',k'), 

Bt =1 (k,k', |k-k'|) -+ - Cl Ap c ( fc )P c ( fc ') 

Bf= 2 (fc,fc',|k-k'|) ->■ 2c 2 P 2 (k-k!)P c {k)P c (k'). (8) 


Note that for the L = 1 case, taking the squeezed limit, the explicit angle dependence V\ (k • k') vanishes because of the property 
of odd function, and only aO(^) term remains. Taking into account the squeezed limit of the other bispectra considered in this 
paper, the large-scale behavior of the transfer function, oc k 2 , and assuming scale invariance of Pq and B £, we analytically 
estimate the following scaling relations on large scales as: 5 


f k~ 2 

k° 


A b oc < 


k - 1 

k - 1 

k~ x 


: local or Co 
: equilateral 
: orthogonal 
: folded 

: ci 


(9) 


Figure [T] shows that our numerical results match these analytic expectations (on large scales). Such estimation also yields 
Ab 0 for the C 2 case, because an angular integral in the convolution 0 is reduced to i n the squeezed limit and 

this becomes zero. This essentially explains the smallness of our C 2 result described in Figure [I] The large-scale amplitude of 
A b in the angle-dependent 0 ^= 0 , 1,2 models decreases as the Fegendre multipole L increases, while, interestingly, they are close 


5 Notice that for the equilateral, orthogonal and folded shapes we are using the standard template shapes 0167). In fact these templates do not recover the 
precise squeezed limit of the theoretical models on which they are based, a limit to which the scale dependent bias can be particularly sensitive to (see, e.g. 
discussions in |68i 69] Despite of this we have used these templates, both because our main goal is to focus on the angle-dependent non-Gaussianity, and 
also to easily compare with previous results in the literature (the use of the standard templates can be anyway taken as examples of the efficiency to use the 
scale dependent halo bias to constrain these shapes). On the contrary, for the scale dependent non-Gaussianity we are using the bispectra predicted by the 
underlying physical model(s), including the large-scale scaling k~ x for the ci shape in the squeezed limit. 
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FIG. 1: Scale-dependent correction to the bias, for different shapes of non-Gaussianity at z — 1. Here we fix the non-linearity parameters as 


.ploc _ f 

J NL — J] 


— /nl — /nl = ci = C2 = 1, and adopt M — 10 11 M o (top left panel), M — 10 12 M© (top right panel), M — 10 13 


M© 


(bottom left panel) and M — 10 M© (bottom right panel). 


to each other on very small scales, due to the difference of their scaling as A b L=0 oc k~ 2 , A b L=1 oc and A b L=2 oc k°. 
We will investigate the observational consequences in the rest of the paper, and in particular in Section VI Our results for the 
standard four shapes are also in excellent agreement with previous studies ED ED, confirming the validity of our numerical 
computations. 

One can see the redshift dependence of A b for each non-Gaussianity model in Figure [2] It is confirmed there that at higher 
z, | Aft| grows while maintaining the overall shapes, because both A\ and A 2 are increasing functions with redshift, on the 
large mass scales that give the best constraints (M ~ 10 14 M o ) [61]. This fact is expected to improve the sensitivity to the 
non-Gaussianity parameters in the analysis with high-redshift data. In particular, we can see that e.g. in the equilateral case, the 
scale-dependent contribution grows with k and with z, so deep surveys that will include both high-z and high-k (because of the 
larger k for which the density field starts to be non-linear at higher redshift) measurements will sensibly improve constraints on 
this shape (see also (60]). 


III. GALAXY SURVEYS 

We investigate constraints on various shapes of non-Gaussianity that can be set by future large-scale galaxy surveys. In 
this paper we are interested in a proof-of-principle investigation to understand if future LSS surveys will allow competitive 
measurements of different shapes of primordial non-Gaussianity. 

Measurements of non-Gaussianity have been done (and forecasted) for spectroscopic (see e.g. 151181145115511701172)) and radio 
continuum surveys (see e.g. [4Jj ( 42, 47 , 73)). In the first case the advantage is the high precision of the redshift information, 
while in the latter it is easier to survey wide areas of the sky and reach higher redshift. Given that constraints on non-Gaussianity 
are in most cases better for large halo mass, scale and redshift probed, and that in general the effect of non-Gaussianity is larger 
at smaller k , so it is not necessary to have very precise redshift information [60], we will focus on forthcoming radio surveys and 
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FIG. 2: Scale-dependent correction to the bias, for different shapes of non-Gaussianity. The relevant input parameters are fixed to the values 
in Figure [l] The results plotted here are for M — 10 14 M 0 . 


assume that a combination of different surveys in optical and radio, complemented by simulations and modeling, will be able to 
assign redshifts to the objects observed in the planned ASKAP/EMU lf74l and the SKA continuum lf49l surveys. In both cases 
we assume a survey observing 30,000 deg. 2 , with total number of objects of ~ 7 x 10 7 and ~ 3 x 10 7 for EMU (flux limit of 
50 yJy, 100 fiJy , respectively), x 10 8 and x 10 9 for SKA with a flux limit of 5 fiJy and 100 nJy , respectively. 

In particular, clustering information can be used to infer redshift distributions. Methods using cross-correlations of samples 
with unknown redshift distributions against photometric and spectroscopic datasets have been proposed [75 , 76 ] and applied to 
real data G7H79), including the FIRST radio survey [80]] . This technique was used also to measure the redshift distribution of 
resolved far infrared HerMES El and lf82l cross-correlated Planck maps against quasars from SDSS to estimate the intensity 
redshift distribution of the Cosmic Infrared Background. A study of the impact of adding redshift information to radio continuum 
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surveys can be found in (83l . 

For the EMU and SKA continuum surveys we use the predicted redshift distributions, bias and halo masses of I49l l73l[84ll86l . 
In particular, we assume that we will be able to distinguish different types of objects, and we follow [ 861 and assign halo masses 
as: 

• Star forming galaxies (SFR): M = 10 11 h~ 1 M Q 

• Starbursts (SBG): M = 5 x lO^/i^M© 

• Radio Quiet Quasars (RQQ): M = 3 x 10 12 /i _1 M q 

• Radio loud AGN (FRI): M = 10 13 

• Radio loud AGN (FRII): M = 




FIG. 3: Redshift distribution for the ASKAP/EMU survey and bias for the different populations observed by radio surveys. 
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FIG. 4: Redshift distributions for the different populations observed by the SKA survey with detection threshold S=5/iJy (left panel) and 
lOOnJy (right panel). 


More details on the different cases studied are in Section [V] 
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IV. POWER SPECTRUM 


We model the galaxy power spectrum as (see e.g. ED): 


Pi(k,»,z) 


[ggM! EM 

[D a (z)} 2 H rei (z) 


[&NG(fc, z) + f(z)/z 2 ] 2 P r 5 {k, z) Fo G(k, n, z) S(k, n, z ) + P shot (z ), 


( 10 ) 


where the superscripts r and 5 indicate real and redshift-space, respectively, and the subscripts 5 and g stand for density and 
galaxies respectively; the superscript ref indicates the values of the angular diameter distance and the Hubble parameter in the 
ACDM case. The shot noise contribution is the inverse of the galaxy number density at that redshift, [^(z)] -1 , and the non- 
Gaussian bias is the one defined in Section [ll| The Redshift-Space Distortion (RSD) corrections come from the fact that the 
real-space position of a source in the radial direction is modified by peculiar velocities due to local overdensities; this effect can 
be modeled as 18811891 : 


6 s (k) = (l + fa 2 )6 r (k), 


where (3 = f/b, and the parameter / is defined as the logarithmic derivative of the growth factor: 


/ = 


dlnD 
din a 


( 11 ) 


( 12 ) 


A real data analysis will need to take into account a variety of additional corrections, that will however not affect considerably 
our study. Examples of these corrections include observational calibrations as well as improved theoretical modeling such as 
wide-angle effects [33l ?9Ql4TQ4l . modifications to the observed radial redshift distribution by Redshift-Space Distortions 11051 . 
cosmic magnification HUES), and its binning lfT07l : however all these are second order effects, so we will not include a 
detailed modeling of them. 

It is also worth noting that the effect of a primordial non-Gaussianity of 0(1) would be degenerate with large-scale ef¬ 
fects 1108111101 . so in a real analysis one would need to take into account relativistic effects 1195-199, HOll111 1 1H1 131 . Here we do 
not include them, because modeling them is computationally expensive, and they should not affect our work, as we focus on the 
precision of the measurement and not on the exact value of the non-Gaussianity parameter. A more detailed work on how GR 
and lensing effects can impact the observable effects of primordial non-Gaussianity shapes and if this could affect forecasts will 
be investigated in a future work CHI. For the /^l case, the effect of neglecting relativistic corrections on errors on and its 
running n n Q, see cm We will then assume for the purpose of this work that one can compute these “relativistic corrections” 
and correct for them. 

Given that in this work we investigate primordial non-Gaussianity shapes that do not limit their effects on the very largest 
scales, it is appropriate to model more carefully the moderately large k regime. We then include in our modeling of the observed 
power spectrum the terms FoG, E, to model the so-called Fingers of God CHI and the error in the determination of the redshift 
of sources, respectively. We can write them as: 


FoG = e j P ee (k) dk ; 


(13) 




(1 + z)g z 
H(z) 


(14) 


where cr v is the velocity dispersion, Poo(k) is the velocity power spectrum and cr z is the expected error in the determination of 
z. 


V. RESULTS 

In this Section we present our results. We investigate the forecasted constraints on /nl for different shapes for future radio 
continuum surveys, assuming some knowledge of the redshift distribution of its sources from cross-ID and simulations, and a 
general spectroscopic survey. In the first case, we also assume that we will be able to identify the type of object observed, and 
this will allow the implementation of the so-called multi-tracer technique 1 1171 , 1181 . Given that we want to focus on future 
constraints on anisotropic models, we also study how what configuration of future surveys could improve constraints on those 
models, and in particular what will be required to reach cr(/^ =1 ) ~ 10 or better, which is comparable to typical values predicted 
by the models l20l[261. 
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Given the specifications of a survey, the Fisher matrix analysis allows us to estimate the errors on the cosmological parameters 
around the fiducial values (see e.g. H19Ul2lh . We write the Fisher Matrix for the power spectrum in the following way: 


P % max pk mg « P + 1 

F a p = / dz dk dji 

** ^min k rn i n •' 1 


V e fi(k, /i , z)k 2 dP{k , /i, 2 :) dP(k, fi, z) 


B n \. 


(15) 


Sir 2 [P(k,^z)} 2 dd a ddp 

where d a ((3) is the a(/3 )~th cosmological parameter, and the effective volume of the survey in the z-th redshift bin is defined as: 

n g {z)P{k,n,z ) 12 


V eff (k, ( 1 5 2 ^) 


_1 +72p(2:)P(fc,/i,2?) 


( 16 ) 


14 is the volume of the survey and n g is the mean comoving number density of galaxies. The last term in Eq. (15) accounts for 
non-linearities induced by the BAO peak f87l : 


Pn\ — & 


-/c 2 £i -k 2 


2 (Ef|-E 2 , 


(17) 


and E_l = E 0 I}, E|| = E 0 (l + f)D , where E 0 is a constant phenomenologically describing the nonlinear diffusion of the BAO 
peak due to nonlinear evolution. From N-body simulations its numerical value is 12.4 h -1 Mpc and depends weakly on k and 
cosmological parameters 1 1221 1231. 

Eq.@ involves an integral over the wavenumber fc; for the maximum scale used it depends on the redshift bin and on the 
geometry of the survey, k m i n = 27rU -1 / 3 , while for & max , following f72l we use the value that gives a value of the variance 
a 2 = 0.36. To be completely correct, the value of k m i n should be computed in a more complicated way that includes detailed 
considerations on the survey geometry, the radial depth of the redshift bin and large-scale modeling (see (EH); however, we 
find that a more complicated modeling of it does not considerably affects our results. We will in any case investigate below the 
effect of varying the choices for the k range used on the measurements of /nl for the anisotropic models. 

We present our results for a Fisher analysis focusing only on the /nl parameter, assuming knowledge of cosmological param¬ 
eters describing e.g. dark energy and curvature, and we allow an overall amplitude of the gaussian bias bo (that might include 
other amplitude uncertainties) to vary. Real constraints might be affected by the uncertainties on other cosmological parameters, 
but here we focus on errors on the scale-dependent bias, and we argue that future priors from CMB and other LSS experiments 
will mitigate the effect on larger Fisher matrices. 

We make use of the so-called multi-tracer technique, that allows the reduction of cosmic variance errors for surveys that 
are observing (and able to discriminate) different types of objects, in particular objects having different bias values. This 
technique was originally proposed in 1117, 1118 1 and then further studied in e.g. 1 125 1: in 1 126 1 it was in particular applied to 
non-Gaussianity analyses. A real data application was performed in fT27l . 

A detailed analysis of how much one can gain using more than one tracer depends on the details of the observed galaxy 
catalog; in particular, the ideal case (see e.g. 11181 ) is when a survey can target a very large number of unbiased sources and a 
subsample of highly biased objects (so a practical way to obtain a more efficient multi-tracer analysis could then be by combining 
a spectroscopic with a photometric catalog). An investigation of the multi-tracer technique for designing a survey measuring 
the local non-Gaussianity parameter /^l was recently presented in l55l . For the practical implementation of this technique we 
follow 1 2 D and we use the 5 populations described in Section [IV| In this case the Fisher matrix becomes: 

f A AT A 7 .dPij(k,fjL,z) r f W1 dPij(k,fJL,z) 

F ap = J dzdkdii 2 ^r(k,/jL,z) -^- [C(k,fi,z) -, (18) 


where T(fc, /i, z) includes the effective volume and the other terms other than the derivatives of Eq.( 15), Pjj = bibjPs + Sij/rii , 
and C is the covariance matrix including auto- and cross- spectra (see £127] for details). 

Here we present our findings when using the radio continuum survey as in Figure [3] modeled after the ASKAP/EMU sur¬ 
vey <741 . for the predicted detection threshold of S = 50/iJy, and for a more conservative case of S = 100/iJy. We then investigate 
the improvements in the constraining power when using an SKA-like survey, assuming a detection threshold of S = 5/iJy, and 
an optimistic S=100 nJy (see |49) for details). A P(k) analysis using radio continuum data will require cross-identification of 
sources in order to assign a redshift to the objects observed. This can be possible by the combination of data with optical surveys 
and the modeling of the redshift distribution of sources using simulations. We compute constraints on the different shapes of 
non-Gaussianity when assuming an EMU survey such as the one of [73] 85], and the two SKA cases, considering three different 
assumptions for the redshift information that will be available for the radio sources detected. We consider: 


a Conservative case, where we divide the catalog into 3 redshift bins, 0<z<l,l<z<2, 2 < z < 6, with cr(z) =0.1; 


• a Optimistic case, where we divide the catalog into 5 redshift bins, 0<z<0.5,0.5<2<l,l<z<2,2<2:<3, 
3 < z < 6 , with cr(z) = 0.03; 
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• a Futuristic case, where we divide the catalog in bins of Az = 0.5 in the range 0 < z < 6 , with a(z) = 0.001. 

In Table [I] we show our results for the above cases. It can be seen that the quality of the redshift information available makes a 
huge difference. 

Errors on the vector fields coefficients will be competitive with CMB ones only for the optimistic SKA 5/iJy case, or a survey 
comparable to SKA 100 nJy will be needed. For the other shapes, we can see that galaxy surveys can more easily be competitive. 

In particular, the EMU survey for the conservative case would not be competitive in constraining non-Gaussianity shapes; in 
the optimistic case, on the other hand, constraints would be comparable to constraints coming from the CMB for the 100/iJy 
survey. The 50/iJy case, instead, would already represent an improvement in measuring /nl- Having a result comparable to the 
CMB one would already be a very good result, as different measurements on different wavelengths, redshift, and instruments 
would represent an important check of Planck results. It is interesting to note that, in the 5 bins case, going from S=100 to 50 
fi Jy would allow the EMU survey to surpass Planck in the precision of /nl measurements. 

For the SKA surveys, even in the conservative case with the 5/iJy detection limit, results are already very interesting and show 
how LSS can realistically improve upon CMB constraints in the not so far future. In the more advanced case of the SKA with 
a 100 nJy detection limit, predicted constraints are impressive and they would allow very stringent tests of inflationary models. 
The futuristic case analyzed here is not necessarily related to radio continuum surveys, but it is a sort of lower limit for the errors 
that this type of analysis could obtain with a future survey. In this case, we model the survey specifications after the SKA survey 
with S=100 nJy, but with redshift information precise enough to divide the catalog in bins of Az = 0.5. It represent more of 
a proof-of-principle investigation, and results are intended to show that in principle LSS can allow extremely precise tests of 
non-Gaussianity of any shape. We also note that our predictions, considering different approaches in the analyses, are in good 
agreement with predictions for the /^l model of e.g. 02[ 47.(861 112 , 1283 . 


Shape 

EMUj, 00 ^ EMUj, 00 ^ 

EMU^ 0,ij!/ 

EMUo 0,ij!/ 

SKA^ SKA^ 

SKA^ 3 ' 

SKAo 00 " j!/ 

Futuristic 

CMB 

ct(/nl) local 

11.94 

5.54 

9.26 

4.37 

1.62 

1.06 

0.67 

0.51 

0.21 

5.7 

ct(/nl) equilateral 

221.14 

79.84 

179.03 

62.58 

22.24 

9.09 

6.18 

2.83 

0.42 

70 

ct(/nl) orthogonal 

102.97 

39.04 

82.25 

30.69 

15.30 

7.40 

6.71 

3.35 

0.54 

33 

ct(/nl) folded 

151.48 

56.45 

121.50 

44.35 

20.29 

9.25 

8.15 

4.14 

0.96 

65 

cr(c L =i) 

1916.29 

721.15 

1519.8 

558.35 

200.62 

78.32 

41.81 

17.54 

2.14 

103 

cr(c L = 2 ) 

10874.9 

4113.82 

8436.7 

2952.5 

1098.93 

393.60 

193.48 

76.55 

9.22 

26 


TABLE I: Results for predicted EMU (with detection threshold S=100/xJy and 50/Uy) and SKA-like surveys (with S=5/iJy and 100 nJy), 
assuming the distribution of sources and their bias as in Figures [5] [4] Results are for the conservative (C) and optimistic (O) cases. We 
also compare with a futuristic survey, representing an upper limit for constraints from this type of analysis. See text for more details. For 
comparison we here depict the CMB temperature limits (68%CL) obtained by the Planck collaboration 0. 


VI. FUTURE CONSTRAINTS ON VECTOR FIELDS MODELS 

Given that the scope of this paper is mainly to investigate if LSS can set meaningful constraints on angle-dependent non- 
Gaussianity, in this Section we try to find the minimum requirements for a galaxy survey to set meaningful constraints on the 
vector fields model parameters cl. In Figures |5j[6] we plot the non-Gaussian correction to the bias for the c\ L= i, 2 parameters 
as a function of redshift, scale and halo mass. It is clear that for the c\ L=1 parameter, very large scales, very high redshifts and 
very large halo masses are required to obtain strong constraints. The c\ L=2 parameter, on the other hand, exhibits a different 
behavior; the amplitude of A b, while decreasing still with halo mass and redshift, decays less dramatically with scale. For this 
reason, as we saw before, increasing the maximum k used in the analysis largely increases the constraining power. 

In the rest of this Section we investigate what are the minimum requirements needed for a galaxy survey in order to reach a 
precision on measurements of c\ L=1 « 10 , which can be considered as a typical (in some cases lower) value predicted by the 
scenarios that produce this kind of non-gaussian signature ll20ll26l 6 ; it is clear from the results in Table [T] and Figures [5] [6j that 
it is unlikely LSS will be competitive with CMB experiments in measuring the parameter c\l =2 . 

In Figure [7] we plot the Fisher element information F Cl=iCl=1 as a function of z and the number density n g in redshift bins 
of Az = 0.1, for different halo masses. As we want to understand if and how it will be possible to reach a precision of ~ 10 
in the measurement of cl=i, we need to reach (in the optimistic hypothesis that we will know the bias from a complementary 
measurement), F Cl=iCl=1 > 0.01. It can be seen that redshift bins for 2 < 2 don’t contribute with a significant amount of 


6 Notice that in the case of the inflationary models proposed in (26) the natural values predicted for ci are typically larger than 10. 
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FIG. 5: Non-Gaussian correction to the bias as a function of scale k and halo mass, for the C\l=i (left panels) and c\ L=2 (right panels) 
parameters. Top panels are for z = 2, bottom panels for z = 7. 


information, so one could use at most ~ 20 bins, and so each F needs to be of the order of « 5 x 10 -4 . In the Figure we plot 
contours of Log(F CL=lCL=1 ) = —4.5, —4, so we indicate what minimum values of n g are needed in order to have the required 
precision. 

The plots show that it will be difficult to reach the precision we wanted in measuring this parameter with a wide spectroscopic 
survey, even at high redshift. It is apparent how a very large number of sources at high redshift would be required in order to 
reach our goal of cf(cl= i) = 10. We then investigate a few of examples of the multi-tracer technique and see how much it can 
be gained with that. We consider the cases of distinguishing between two populations with mass Mi < M 2 , and assume that n 2 
is an order of magnitude smaller for each order of magnitude that M increases. We show our results in Figure [8] and we can 
see that the multi-tracer case helps considerably. It is still required to have a survey with a high number density at high-z, but 
the requirements are obviously less stringent, and potentially reachable. In the plots we show, as in the single tracer case, the 
Fisher element information F Cl=lCL=1 as a function of z and the number density n g9 for two different cases of multi-tracer, for 
a wide (30,000 deg. 2 ) and a narrower (5,000 deg. 2 ) surveys. We overplot lines showing the minimum number density required 
for the smaller halo mass population, per redshift bin of Az = 01., in order to reach ct(cl= i) ~ 10. We can see that for this 
parameter f s k y is not as important as it is for the more studied local shape /nl (see e.g. 145)), and so a deep survey targeting a 
large number of highly massive objects would be the easiest way to measure cl =i with high precision. 

It is beyond the scope of this paper to investigate details of future surveys and how to reach the required numbers; we leave 
this as a guide for designing future experiments and a future study. 


VII. CONCLUSIONS 

In this paper we investigated the possibility of using large-scale radio surveys to measure primordial non-Gaussianity in 
configurations other than the widely investigated “local” shape. In particular, we focused on the possibility of constraining 
vector field models with the large-scale structure of the Universe. This has not been done so far, and the advent of large galaxy 
surveys in the next few years makes this timely. 
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FIG. 6: Non-Gaussian correction to the bias as a function of redshift z and halo mass, for the c\ L = 1 (left panels) and c\ L=2 (right panels) 
parameters. Top panels are for k = 0.001, bottom panels for k = 0.1. 


We forecasted the precision on measurements of the primordial non-Gaussanity parameter /nl for the local, equilateral, 
orthogonal and folded shapes, and the vector models described by the c\ L=1 ^ coefficients, coming from the galaxy power 
spectrum measured by future galaxy surveys. We computed constraints that will be possible to obtain by having some redshift 
information added to forthcoming radio continuum surveys such as EMU and the SKA, in different configurations. We found that 
an optimistic but realistic analysis could allow measurements of the /nl parameter competitive or improving the current limits 
coming from CMB experiments. In particular, we found that having the ability to divide the galaxy distribution function into 5 
bins, already the coming EMU survey could give constraints better than current Planck limits. If we assume that the distribution 
can be divided only into 3 bins, the SKA will be needed in order to improve upon CMB measurements. In a more futuristic 
perspective, as a proof of principle investigation, we found that LSS analyses could be extremely powerful in measuring non- 
Gaussianity parameters from any shape, provided that redshift information for a large number of galaxies over a large fraction of 
the sky and at high-z will be available. Finally, we studied the specifications required for a generic future galaxy survey in order 
to reach a precision of ct(cl= i) = 10, that would allow test of the primordial Universe models with parity-violating anisotropic 
sources, e.g., a U(l) vector field coupled to a pseudoscalar and large-scale helical magnetic fields in radiation dominated era. To 
achieve an even better precision on measurements of non-Gaussianity parameters, it will be necessary to have a high-redshift 
wide survey and the use of higher order correlations such as the bispectrum (see e.g. E3133133). 
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